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Optimal Neumann control for the ID wave equation: Finite 
horizon, inhnite horizon, boundary tracking terms and the 

turnpike property 

Martin Gugat* Emmanuel Trelat^ Enrique Zuazua^ 


Abstract 

We consider a vibrating string that is fixed at one end with Neumann control action at 
the other end. We investigate the optimal control problem of steering this system from given 
initial data to rest, in time T, by minimizing an objective functional that is the convex sum 
of the L^-norm of the control and of a boundary Neumann tracking term. 

We provide an explicit solution of this optimal control problem, showing that if the weight 
of the tracking term is positive, then the optimal control action is concentrated at the beginning 
and at the end of the time interval, and in-between it decays exponentially. We show that the 
optimal control can actually be written in that case as the sum of an exponentially decaying 
term and of an exponentially increasing term. This implies that, if the time T is large the 
optimal trajectory approximately consists of three arcs, where the first and the third short-time 
arcs are transient arcs, and in the middle arc the optimal control and the corresponding state 
are exponentially close to 0. This is an example for a turnpike phenomenon for a problem 
of optimal boundary control. If T = -|-oo (infinite horizon time problem), then only the 
exponentially decaying component of the control remains, and the norms of the optimal control 
action and of the optimal state decay exponentially in time. In contrast to this situation if the 
weight of the tracking term is zero and only the control cost is minimized, then the optimal 
control is distributed uniformly along the whole interval [0, T] and coincides with the control 
given by the Hilbert Uniqueness Method. 

Keywords: Vibrating string, Neumann boundary control, turnpike phenomenon, exponential 
stability, energy decay, exact control, infinite horizon optimal control, similarity theorem, receeding 
horizon. 


1 Introduction 

The turnpike property has been discussed recently for the optimal control of linear systems gov¬ 
erned by partial differential equations, see [14]. Turnpike theory has originally been discussed in 
economics, see [16]. For systems governed by ordinary differential equations the turnpike property 
has been discussed for example in [17], also for the nonlinear case. The turnpike property states 
loosely speaking that if the objective function penalizes both control cost and the difference of the 
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optimal trajectory to a given desired stationary state, the optimal controls will steer the system 
quickly to the desired stationary state and then the system will remain on this path most of the 
time. In Section 4 of [14], optimal control problems with the wave equation are considered on a 
given finite time interval [0,T]. The control is distributed in the interior of the domain and no 
conditions for the terminal state at the time T are prescribed. In this paper we consider problems 
of optimal boundary control of the wave equation. We consider both problems with infinite time 
horizon and problems with finite time horizon. For the finite time problems, we prescribe exact 
terminal conditions. If the objective function only penalizes the control cost, for the ID case in 
[6] it has been shown that the optimal controls are periodic, so in particular they do not have the 
turnpike structure. This illustrates that the turnpike property heavily depends on the choice of 
the objective functions, that has to couple the control cost with the penalization of the distance 
to the desired state. 

We consider a system governed by the one-dimensional wave equation on a finite space interval, 
with a homogeneous Dirichlet boundary condition at one side, and a Neumann boundary control 
action at the other: 

dttyit,x) = d^^y{t,x), (t,x) G M X (0,1), 

y{t, 0) = 0, d^yit, 1) = u{t), t G K, 

where the control u belongs to the class of square-integrable functions. 

Let yo G 1) be such that yo(0) = 0, and let yi G L^(0,1) be arbitrary. For every 

u G L^(0, -boo), there exists a unique solution y G -boo; H^{0, 1)) nC'^(0, -boo; L^(0,1)) of (1) 
such that y(0, •) = yo{-) and dty{0, •) = yi{-). 

As is well known, the system is exactly controllable if and only if T ^ 2 (see for example [7]). 
In this paper, given any A G [0,1] and any T G [0, -boo], we consider the optimal control problem 
(OCP))]’ of finding a control u G L^{0,T) minimizing the objective functional 

J\{u)= j {{I - X) {d^y{t,Q)f‘^ dt, (2) 

such that the corresponding solution of (1), with y(0,-) = yo(-) and dty{0,-) = yi{‘), satisfies 
y{T, •) = dty{T, •) = 0 at the final time (exact null controllability problem). If T = -boo, then one 
can drop the final constraint requirement, which, by the way, happens to be automatically satisfied 
in the sense of a limit (this is a well known result coming from Riccati theory). 

For A = 1, this very classical optimal control problem (minimization of the norm of the 
control) has been considered, e.g., in [10, 12, 15] (see also [2] for optimal control problems consisting 
of satisfying consumer demands at the boundary of the system, such as in gas transportation 
networks), and it is easy to see that the optimal control is periodic (see [6]), with a period equal 
to 4, that is, twice the time needed by a wave starting at the boundary point where the control 
acts to return to that point. Note that, in this case, the optimal control is as well given by the 
Hilbert Uniqueness Method (see [12]). 

If A < 1, then the objective functional involves a nontrivial boundary tracking term. This 
tracking term may be considered as a boundary observation of the space derivative of the state 
at the uncontrolled end of the string. As we are going to prove, in that case, the optimal control 
action is then essentially concentrated at the beginning and at the end of the time interval [0,r]. 
More precisely, the optimal control can be written as the sum of an exponentially decaying term 
and of an exponentially increasing term. 

As a consequence, if T is large then the optimal control, solution of (OCP))^’, approximately 
consists of three pieces: the first and the third pieces are in short-time, and are transient arcs; 
the middle arc is in long time, and is exponentially close to 0. This is a turnpike phenomenon 
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(see [14, 17]), meaning that the optimal trajectory, starting from given initial data, very quickly 
approximately reaches the steady-state ( 0 , 0 ) (within exponentially short time, say e > 0 ), then 
remains exponentially close to that steady-state within long time (say, over the time interval 
[e, T — ej), and, in the last short-time part [T — £, T], leaves this neighborhood in order to quickly 
reach its target. 

In this approximate picture, if T = -|-oo (infinite horizon), then the last transient arc does 
not exist since the infinite-horizon target is the steady-state (0,0). In that case, the norm of the 
optimal control decays exponentially in time, and the same is true for the optimal state. Indeed, 
smallness of the observation term for a sufficiently long time interval with zero control implies 
proportional smallness of the state (this follows from an observability inequality, see [ 20 ]). 

Another possible picture illustrating the turnpike behavior is the following. For T large, the 
optimal trajectory of (OCP)J approximately consists of three arcs: the first arc is the solution 
of (OCP)^ (infinite horizon problem), forward in time, and converges exponentially to 0. The 
second arc, occupying the main (middle) part of the time interval, is the steady-state 0. The third 
arc is the solution of (OCP)^, but backward in time. Note that the optimal control problem 
(OCP)^ fits into the well known Linear Quadratic Riccati theory. 

In all cases, we will provide completely explicit formulas for the optimal controls, which explain 
and imply the turnpike behavior observed for A < I. This is in contrast with the case A = 0 for 
which the control action is distributed uniformly along the time interval [0, T], In addition, we will 
also establish a similarity theorem showing that, for every T that is a positive even integer, there 
exists an appropriate weight A < I for which the optimal solutions of (OCP))]’ and of (OCP)“ 
coincide along [ 0 , 2 ]. 

In this paper we focus on problems governed by the ID wave equation. In order to illustrate 
the generality of the turnpike phenomenon, before we turn to the ID wave equation we consider 
an example in a more general framework with a strongly continuous semigroup in the following 
section. 


2 A general remark 

Let a Hilbert space X be given. Let A be the generator of a strongly continuous semigroup {Tt)t^o 
(for the definitions see for example [11], [18]). Let U be another Hilbert space that contains the 
controls. Let the linear operator B : U ^ X he given. Let a time T > 0, a weight 7 > 0 for the 
control cost and an initial state yo & X he given. For u S L^((0,T); U) and t S [0,r], define 

'&((«)= [ Tt-aBu{a)da. 

Jo 

We consider a system that is governed by the differential equation y'{t) = A y{t) B u{t), t S [0, T] 
with the initial condition y{0) = yQ £ X. We assume that there exists a time Tmin > 0 such that 
for all T ^ Tmin the considered system is exactly controllable in the sense that 


Ran = X. 


Let us assume that T > T^in, then there exists a control function u £ L^((0,T); U) such that the 
terminal constrain y{T) = 0 holds. Let A G (0,1) be given. Consider the problem of optimal exact 
control 

f min/J’(l - X){y{t),y{t))x + \{u{t),u(t))u dt 

■\ subject to y'{t) = Ay{t) + Bu{t),t £ ]^,T]] (3) 

I y( 0 ) = 2/0, y{T) = 0. 
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The static optimal control problem corresponding to (3) is 


f min(l - A)(j/, 2 /)x + A(u,m)c/ 
[ subject to 0 = Ay + Bu. 


Obviously the solution of (4) is zero. The solution of the static optimal control problem (4) 
determines the turnpike which in our case is (utp,ytp) = (0,0). Results about the convergence of 
the long time average of the optimal controls to the turnpike control can be found in [1], [19]. 

In order to determine the structure of the optimal control that solves (3) we look at the 
necessary optimality conditions. For all di S L^((0, T); t/), p S C((0,T);X), 82 G C((0,T);X) 
with ^ 2 ( 0 ) = 0 , 52 {T) = 0 we have 

0 = f X (m(s), di(s ));7 + (1 - X){y{s), 62 {s))x + - ^i 52 (s) - BSi{s), p{s))x ds. 

Jo 

This yields the optimality system 

( y' = Ay + Bu, 

{ p' = -A*p + y, ( 5 ) 

[ U = J^B*p 

with the conditions y{0) = yo, y{T) = 0. Hence we get 

— Ap' = AA*p — Ay = AA*p — y' + ^ BB* p. ( 6 ) 

A 

By taking the time derivative in the first order equation for p we get the second order equation 
p" = —A*p' + y'. This implies 

y'=p" + A*p'. (7) 

Now we can use (7) to eliminate y' from ( 6 ) and get 

-Ap' = AA*p - p" - A*p' + ^ ^ BB* p. 

A 

This yields 

p" = {AA*+^^BB*)p+{A-A*)p'. (8) 

A 


2.1 Skew-adjoint operators 

Now we consider the case where A is skew-adjoint, that is A* = — A. In [14], turnpike inequalities 
for the case of the wave equation where A* = —A are given in Section 4. Equation ( 8 ) yields 

p" = Lp + 2Ap' (9) 


where L = AA* + ^ BB*. 

If there exist solutions iA_|_, of the operator equation 

D'^ =L + 2AD 


this yields solutions of the form 


p{t) = exp(iA+ 1) Pi + exp(iA_ t)p 2 
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where Pi, P 2 & X are chosen such that for the state y = p' + A*p we have y(0) = 0 and y{T) = 0. 
Note that L + is positive in the sense that (x, (L + A^)x)x ^ 0 for all x G X. If 

ABB* is skew adjoint (for example if B is the identity) we have D± = A± (L + . 

Now we assume that A and B B* are diagonalizable with the same sequence of orthonormal 
eigenfunctions {‘Pk)k and that the real parts of the eigenvalues of B B* are bounded from below 
by > 0. Then (9) yields a sequence of ordinary differential equations for hk{t) = {p(t),(pk)x 
namely 

h'l, = {ifk, Lpk)x hk + Aipk)x h'k- (10) 

With the roots <5^, of the characteristic polynomial 

Pk{z) = - 2{(pk, A(pk)x z - {ipk, Lipk)x 


we get the solutions 

hk{t) = ■Ufeexp((5^t) + Vkex.p{S'^t) 

where Re(^^) ^ w > 0 and Re(<5^) < —oj < 0. The coefficients Uk and Vk are chosen such that 
p'(0) = 2/0 + Ap{0) and p'{T) = ApiT), because then we have y(0) = //o and y{T) = 0. In fact, 
this yields a constant Cmin that is independent of //o and T such that for all T ^ Tmin, we have 
the inequality 

X] + |exp(2(5^T)| lufcl^) ^ Cmin II2/0IIX- 

k 

Using Parseval’s equation we get the inequality 


^ e 


+e |exp(2(5+T)| lufcp)^ 

< + 2v^||2/o||x. 


1/2 


( 11 ) 

( 12 ) 


Inequality (11)-(12) is a turnpike inequality for p. It states that the norm of p(t) is bounded 
above by a sum of a part that is exponentially decreasing with time and a second part that is 
exponentially increasing towards T. The optimal control has the form u = ^-^B* p, so it also 
shows a turnpike structure. Note that the optimal control norms are decreasing with T, hence 
they are uniformly bounded. 


2.2 Self-adjoint operators 

Now we consider the case that A is self-adjoint. In [14], turnpike inequalities for the parabolic case 
where A* = A are given in Section 3. Equation (8) yields 

p" = Lp (13) 

where L = AA* + ^BB*. Thus we get 

1 11 
p(t) = cosh(tL2)p(0)-I-L 2 sinh(tL2)p'(0) 

with the cosh and sinh operators as defined in [4]. For the optimal state we have 

y = p' + A*p. 

The equations 2/(0) = yo, y{T) = 0 yield a system of linear equations for p(0), p'(0). In particular 
we have p'(0) = yo — A*p{0). 
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Now let us assume that L is diagonalizable and the eigenvalues of L are bounded from below 
by uP' > 0. Then we have the turnpike inequality 

IbWII < ^exp(-wt) (IIp(O)II + |lL-^/V(0)||) + ^exp(-w(r-t)) 

(II exp(Li/2T)p(0)|l + II exp(Li/2r)L-i/V(0)|l). 


The optimal control has the form 


u(t) 


-B* cosh(t L2 )p(0) + —B* 
7 7 


1 1 

'2 sinh(tL2y(0). 


(14) 


This means that also the optimal control can be represented as the sum of families of increasing 
and decreasing exponentials with rates lo. 


3 The main results 

Now we come to our results about optimal control problems for a system governed by (1). Let 
j/o G -ff^(0,1) be such that yo(0) = 0, and let yi S L^(0,1) be arbitrary. 


3.1 Explicit formulas 

We have the following result, giving the explicit solution of (OCP)T for any A G [0,1] and any 
T G [2, +oo]. Note that, when T = +oo, we have to assume that A < 1 to ensure well-posedness. 


Theorem 1. For every T G [2,+oo] and every A G [0,1], the problem (OCP))^ has a unique 
optimal control solution denoted by rtj. 

1. We assume that T = 2n for some n G N*. 

• // A = 1, then the optimal control uf, solution of (OCP)J' is 2-anti-periodic, and thus 
Aperiodic, meaning that 

uj{t-\-2k) = {-!)'"Ui{t), (15) 

for every t G (0, 2) and for every fc G N such that t -\- 2k ^ T, and moreover, 


(t) = 


^ - t) - yii^ - t)), tG{0,l), 

^ (2/ob-1)+ ?/ib-1)) > ie(l,2). 


(16) 




7/ A < 1, then the optimal control solution of (OCP) J is the sum of an exponentially 
decaying term and of an exponentially increasing one. More precisely, defining the real 
number z\ G (—1,0) by 


-A 

2-A + 2v/T^’ 


(17) 


we have 

ul{t + 2k) = zlf+{t) + \f-{t), 


for every t G (0, 2) and every fc G N such that t + 2fc ^ T, where 


(18) 


hit) 


1 -\- Z\ 


F{t-1), 


f_{t) = l^Fit-l), 

J- - 9 r) 


(19) 
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with F G L^(—l,l) defined by 


m = 


5 ( 2 / 0 W +2/1W) ) 


tG [0,1). 


(20) 


2. We assume that T = +00 and that A < 1. 


If X = 0, then the optimal control u'^, solution o/(OCP)-, coincides along the time interval 
[0,2] with the optimal control u\, solution o/(OCP)i. 

If 0 < X < 1, the optimal control solution of (OCP)“, is given along the time interval 
[0,2] by 


uTit) = 

and moreover, we have 


l+^A (,,! 

2 ^ (2/0 


( 2 / 0 ( 1 -i)-2/i(l-i)), 1G(0,1), 


( 2 / 0 ( 1 -1)+ 2/1 (1-1)) , 1G(1,2), 


( 21 ) 

( 22 ) 


u^{t + 2k) = zlufit), 
for every t G (0, 2) and every k gN* . 

The corresponding optimal state yff decays exponentially, in the sense that there exists Cq > 0 
such that 


(^{9xyTit + 2k, x)f + {dtyT{t + 2k, x)f'j dx 


^ Cq\zx 


\2k 


+yii^)'^) dx, (23) 


for every t G (0, 2) and every A: G N*. 


Theorem 1 is proved in Section 4.2. 

For A = 1, that is, when there is no tracking term in the objective functional, the explicit 
solution of (OCP)^ given above has already been computed in [6, Theorem 2.1]. In this case, the 
problem consists of minimizing the norm of the (Neumann) control. The optimal control uj, 
whose explicit formula is given above, can also be characterized as well with the famous Hilbert 
Uniqueness Method (see [12]) and is then often referred to as the HUM control. 

Here, there is no dissipation induced by the objective functional (no tracking term), the optimal 
control is periodic, and is uniformly distributed over the time interval [0, T], in the sense that there 
is no energy decay. 

In contrast, if A < 1, the control is the sum of two terms, one of which is exponentially 
decreasing, and the other being exponentially increasing. For T large enough, this implies the 
turnpike phenomenon, stated in details in Section 3.2. 


Remark 1. For A = 0, the solution of (OCP)[j° coincides with the solution of the problem of 
optimal feedback control studied in [9]. 


Remark 2. The estimate (23) is clearly equivalent to 

({d,,y')f{t,x) f + {dtyT{t,x)f^dx^Cie~^^*J {y^ixf + yi{x)'^) dx, 

for every time t ^ 0, for some positive constants Ci and y not depending on the initial state (see 
also [8, Lemma 2]). 
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Remark 3. It is well known that the solution of the infinite horizon problem (OCP)^ can also be 
expressed in feedback form (Linear Quadratic Riccati theory), see for example [5]. More precisely, 
the velocity feedback 

dxy{t, 1) = I dty{t, 1) 

- 1 

generates the same state as the one generated by the optimal control . 


Remark 4. In the above results, we considered only the null steady-state, but we can easily replace 
it with any other steady-state, as follows. Any steady-state of (1) is given by y{x) = ax, for some 
(T G M. Then, all results therein can be written in terms of such a steady-state: it suffices to 
replace, everywhere, y(t,x) with y(t,x) — ax, and dxy{t,x) with dxy(t,x) — a. For instance, the 
right boundary condition becomes d^yfi,!) = cr + u{t), the final conditions become y{T,x) = ax 
and dty(T, x) = 0, and the objective functional becomes 

Jx{u)= J ({1 - X) {dxyit,0) - af-\-Xu{t)'^^ dt. 

3.2 Consequence: the turnpike behavior 

From Theorem 1 and from the previous discussions, we infer the following consequence on the 
qualitative behavior of the optimal solution. 

Corollary 1. For every X G [0,1), then there exist Ci > 0 and y > 0 such that, for every T ^ 2, 
for all initial conditions {ya,yi) G 11^(0,1) x L^(0,1) with yoiff) = 0, the optimal solution of 
(OCP)J satisfies the estimate 

J (^{dxylit, x)Y-G {dtylit, x))'^'j dx J {yoix)"^-G yi{xf) dx, (24) 

for every t G [0, T]. 

In the estimate (24), what is important to see is that the term is equal to 1 at times 

t = 0 and t = T, but it is exponentially small in the middle of the interval. It becomes even 
smaller and smaller when T is taken larger. This estimate implies the turnpike behavior described 
previously: short-time arcs at the beginning and at the end of the interval are devoted to satisfy 
the terminal constraints, and in-between, the trajectory remains essentially close to rest. 

3.3 Similarity result 

We next state the following similarity result: for any final time T that is a positive even integer, 
there exists a weight A such that the optimal solutions of (OCP)“ and (OCP)“ coincide along 
the subinterval [0,2] of [0,T]. 

Theorem 2. Given any T G 2N*, we choose A > 0 such that 

= I - 1. (25) 


Then we have 


Ui{t) =uff{t), VtG (0,2), 


(26) 



and 


^I’O) ~ ''^?’(’)||i 2 ( 2 fc_ 2 fe+ 2 ) ^ 1^ ~ I^aI^I (ll2/olU^(0,l) + I| 2 /i||l 2 ( 0 . 1 )) • (27) 

for every k £ {0,1,..., (T — 2)/2}. 

The proof of Theorem 2 is done in Section 4.3. 

3.4 Numerical illustration 

We set yo{x) = 4sin(7rx/2) and yi{x) = 0, for every x G [0,1]. From Theorem 1, if 0 < A < 1 then 
the optimal control solution of (OCP))^’ is given by 

ul{t + 2k) = \z\\^ (1 + zx)Tr sin (^^{t + 2fc)^ , 


for t G (0,2) and A: G N. 

The graph of dxyj{t,x) is provided on Figure 1, for T = 20, with A = 24/25 on Figure 1(a) 
and A = 99/100 on Figure 1(b). The control uj){t) = dxyx{t, 1) is the boundary trace at the back. 


8., 
6 . 



0 0 

(a) A = 24/25 and T = 20. 



(b) A = 99/100 and T = 20. 


Figure 1: Plot of dxy{t,x). 


These hgures illustrate that the norm of the optimal state decays faster if A is smaller, as 
expected. However, smaller values of A cause larger oscillations. Note that z;, = —2/3 if A = 
24/25, and za = —9/11 if A = 99/100. Moreover, as pointed out in [6], if T G 2N* then Ui{t + 
2k) = TjiTrsin + 2k)), for all t and k such that t + 2k G [0,T] (see also [6, Figure 4] for the 
corresponding optimal state, with T = 10, up to the factor 27r). 


4 Proofs 

4.1 Well-posedness of the initial-boundary value problem 

Let yo G H^{0, 1) be such that yo(0) = 0, and let yi G L^(0,1) be arbitrary. Let T ^ 2, and let 
u G L^(0, T) be hxed. As a preliminary result, we study the well-posedness of the initial boundary 
value problem (1) for a fixed control u, and with the fixed initial data {yo,yi). The analysis is 
similar to the one done in [8]. 
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We search a solution given as the sum of two traveling waves, i.e., 


y{t, x) = a{x + t) + (3{x - t), 


where the functions a and (3 are to be determined from the initial data and from the boundary 
data. First of all, to match the initial conditions, we must have 

a{t) = ^(^yoit) + yi{s)ds^+Co, (28) 

Pit) = ^(^Voit)-J^yiis)ds^-Co, (29) 

where Co is a real number. Besides, the boundary condition y{t, 0) = 0 implies that 

Pi-s) = -a(s), (30) 


for almost every s > 0. The boundary condition at a; = 1 leads to a'(l + t) = u{t) — /3'(1 — t), and 
integrating in time, we get 

a(t + 1) =/3(1 — t) + f w(s) ds + a(l) —/?(!). 

Jo 

Using (28) and (29), we have a(l) — /?(!) = J^yiis)ds + 2Co, and therefore, choosing Co = 
fo di(s) ds, we get a(l) — /3(1) = 0 and 


+ 1) — P(i ~ t) + 



(31) 


Using (30), the values of a for t S (0,1), given by (28), determine those of /3 for t G (—1,0). The 
values of /3 for t G (0,1) are given by (29). Now, knowing /3 on the interval (—1,1), we deduce 
from (31) the values of a on the interval (1, 3). 

Using (30), we get a(t + 1) = —a(t — 1) + /q u(s) ds, for f ^ 1, or equivalently, 


ait + 2) = —ait) + / w(s) ds, (32) 

Jo 

for t > 0. 

Using (32) enables us to determine a iteratively: starting with a on the interval (1,3), the 
values of u{t) yield those of a on (3,5), and then using (32), we determine a on (7,9), etc. 

In order to express everything in terms of a only (without using j3), we extend the domain of 
a so that it contains (—1,0). We get the values of a on (—1, 0) by using (30) for s G (0,1), which 
yields ait) = —/3(—t) for t G (—1,0) with the values of p on (0,1) given by (29). Then, using (31), 
we get ait + 2) = —ait) + 'u(s) ds for t G (—1,0). We have the following lemma. 

Lemma 1. Let yg G 77^(0,1) be such that yoiO) = 0, and let yi G L^(0,1) be arbitrary. We set 


and we define a G L^(—1,1) by 


Co = -- 


yiis)ds. 


(33) 


ait) = 


k {-yoi-t) + Jq * yiis) ds'^ + Co if fG(-l,0), 
5 (do (7) + fo 2/1 (s) ds'j +Co if t e [0,1). 


(34) 
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Let T = 2K G N*, and let u G i^(0, T) be fixed. 

The function a, defined by iteration according to 

a{t + 2k) = —a{t + 2{k — 1)) + / u{s) ds, (35) 

Jo 

for every t G (—1,1) and every k G N such that t < T + I — 2k, is well defined on the interval 
(—1,T + 1), and belongs to + 1). 

Proof: From the construction, it is clear that a |(fe-i,fc) S H^{k,k + 1), for every k G N. To 
prove that a G T + 1), it suffices to prove that a is continuous. Since a(0+) = a(0“) = Co, 

a is continuous at t = 0. Using (34), a is continuous as well on (—1,1). 

At t = 1, using (33) we get a(l+) = 0(1”) = |yo(l), and hence a is continuous at t = 1. Then, 
at this step, we have obtained that a is continuous on (—1,3). 

We then proceed by induction. Let k G N*. We assume that a is continuous on the interval 
(—1,1 + 2k). Then a((—1 + 2fc)“) = Q!((—1 + 2^)+). Using (35), we have 

p2k p2k 

a((l + 2fc)—) = —a((l + 2(fc — 1))“) + / u(s) ds = —a((—1 + 2A:)“) + / u{s) ds 

Jo Jo 

p2k 

= —a((—1 + 2fc)''") + / ii(s) ds = a((—1 + 2(fc + 1))"*") = q ;((1 + 2^)“*"). 

Jo 

Since a is defined by (34), we infer that a is continuous on (—1, l + 2(fc + l)) for fc + 1 ^ K. Lemma 
1 is proved. □ 

Using Lemma 1, we are now in a position to compute the solution of the initial boundary value 
problem under consideration in this subsection. 

Proposition 1. Let j/q € -^^(0,1) be sueh that yo(0) = 0, and let yi G L^(0,1) be arbitrary. Let 
T G 2N*, and let u G L^(0,T) be fixed. We consider the function a defined in Lemma 1 by (34). 
Then the solution of (1), associated with the eontrol u and with the initial data (yo:2/i): *5 given 
by 

y{t,x) = a(t + x) — a{t — x), (36) 

for all {t, x) G (0, T) x (0,1). 

Proof: The construction of a implies that y, defined by (36), is a solution of the initial 
boundary value problem under consideration. We conclude by Cauchy uniqueness. □ 

4.2 Proof of Theorem 1 

4.2.1 Case T < +oo 

Let yj be the state generated by the control uj) defined in the theorem. Let us first prove that yj) 
satisfies the terminal constraints 


yf(r,.) = 0, dty^(T,-)=0. (37) 

It suffices to prove that a'(z) = 0, for z S (T — 1, T + 1). From (36), we have yf{t,x) = 
a\(t + x) — a\(t — x), with ax defined by (34). The definition (20) implies that F(t) = a'^(t). 
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Hence, we have 


Ait-^) = YZ 


1 _ 2n 

±nt-i) = 


1-4" 


~2n 


-1 


Fit-1) 


l-z|" 


1- I 1 + Z\ 


4 x/-('). 


where the last equality follows from (19). 

By (35) we have 

+ 1) = ~ 1) + (38) 

for t e (0,2). Using (18), this yields a'^it + 1) = —a'^it — 1) + /+(t) + /-(t), for t € (0,2), and 
then, using (19), 

«a(* + 1) = — f+ii) + , ^ 1 /-(^)- 

1 + zx 1+4 

By induction, thanks to (38) and (18), this implies that 


k — 

rrKi “ 1 + 2fc) = ——/+(t) + —4i-/_(t), (39) 

r + ^A i + 

for every t G (0, 2) and every fc G N such that 2k ^ T. Taking k = T/2, we get 

1 

~n — 

a\it -1 + T) = ^ Uit) + (40) 

^ I -L “r 

Using (19), we infer that 

a4t-i + r)= + = (41) 

and hence the state yj satisfies the terminal conditions (37). 

For a control of the form u = uj) + h, the generated state is y = yj) + yh, where yt is 
the state generated by the perturbation control h, with the boundary conditions yhit,0) = 0, 
dxVhi't^ 1) = li(l)i E^rid null initial conditions. We only consider variations h for which yniT, ■) = 0 
and dtyhiT, •) = 0. Using (36), we have 


yhit, x) = ah{t + x) - Uhit- x). 


Since ?//j(0, •) = dtyhiO,-) = 0, we must have — 1) = ahit) = 0. Moreover, owing to the 
terminal constraints, we must have = 0 along (T — 1,T+ 1). 

The value of the objective functional of (OCP)j(’ is 


((1 - X){dxy{t, 0))^ + Au(t)^) dt 


((1 - A) {ida;ylit,0))'^ + id^yhit,0))‘^ + 2d^ylit,G)d^yh{t,0)) 

+ A + 1^(0^ + 2ma 


dt. 
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We consider the linear part 


Lx{h) = 2 / {{1 - X)dxyl{t,0)dj;yhit,0) + \ul{t)h{t)) dt. 

Jo 

Since 0) = 2a\{t), da;yhit,0) = 2a'^(t), and 

ulit) = d^ylit,!) = a'xit + 1) + a'xit - 1), 

Ht) = d^yh(t, 1) = + 1) + - 1), 

we get 

(8(l-A)aUiKW 

+ 2A {a'x{t + 1) + a'x{t — 1)) + 1) + — 1)) ^ dt 

cl , 

= (8(1 - A)a;,(t ++ j) 

+ 2A {a'xit + 1 + j) + a'xit — 1 + j)) (a(j(t + 1 + j) + — 1 + j)) ^ dt 

T-l ^1 

= 8(1-A)^/ a'^{t + j)a'xit + j) dt 

j=i -^0 
T —2 1 

+ 2A^ / + j) + 2 + j) + + j)) dt 

,=i -lo 

rl 

+ 2A ^ / a'^(t + j) (a^(t + j) + — 2 + j)) dt 

,=i Jo 

T—2 I 

j=i -^0 ^ 

A (®a(^ + 2 + j) + a^;^(t — 2 + j) + 2a^(t + j)) ^ dt 

Defining the characteristic polynomial by 

P\{z) = Xz^ “b (4 — 2A)z + A, (42) 

we have p\{z\) = 0 and px{llzx) = 0. Using (39), we have 

Xa'x{t + 2 + j) + (4 — 2X)a'x{t + j) + Xa'x{t — 2 + j) = 0, 

for every t S (0,1) and every j G {1, 2,..., T — 2}. This implies that Lx{h) = 0. Now, concerning 
the value of the objective functional of (OCP)j^, for any h such that yh{T, 0=0 dtyh{T, •) = 0, 
we infer that 

J ((1 - A) (5a;j/(t,0))^ + Au(t)2^ dt > y ((1 - A) (5a;yf(t,0))^ + AuJ(t)^) dt, 

with a strict inequality whenever h 0. It follows that uj^ is the unique optimal solution of 
(OCP))[', as soon as A > 0. If A = 0, then the result also follows from the representation of Lo(^)- 
However, in this case the characteristic polynomial Po(z) = Az has only one root given by Zq = 0. 
Theorem I is proved for T < +oo. 
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4.2.2 Case T = +oo 

We are going to use the previously established well-posedness results. 

Let y'^ be the state generated by the control uf" defined in the theorem. For a control of the 
form u = + h, the generated state is y = +yh, where yt is the state generated by the control 

h, with null initial conditions and with the boundary conditions yh(t, 0) = 0 and dxyh(t, 1) = h{t). 
The value of the objective functional of (OCP )r is 

/ ((1 - A)(5a:y(t, 0))^ + Au(t)^) dt 

Jo 

= J ((1 - A) {{dxyl{t,0))'^ + {dxyhit,0))'^ + 2dxyT{t,Q)dxyh{t,0)) 

+ A (t)^ + h{ty + 2u^{t)h{t)) ^ dt. 

We consider the linear part 

p+oa 

Lx{h) = 2 / {{1-X)dxyT{t,0)dxyh{t,0) +Xu^{t)h{t))dt. 

Jo 

Using (36) , we have j/“ (t,x) = a\{t + x) — a\(t — x) and yh(t,x) = ah(t + x) — ah(t — x), with a\ 
given by (34), and a/i = 0 on (—1,1). It follows that d^y^^it, 0) = 2a'^(t), d^ytit, 0) = 2a'j^(t), and 

uT{t) = dxVTit, 1 ) = oi'xit + 1 ) + oi'xit - 1 )> 

Ht) = dxVhit, 1) = a'^(t + 1) + - 1), 

and therefore, 

Lx{h) = (8(1 - X)a\{t)a'f,{t) 

+ 2A + 1) + a'x{t — 1)) + 1) + a/j(t — 1)) ^ dt 

+CXD -1 

= (8(1 - A)a';,(t + j)a;,(t+ j) 

+ 2A {a'xit + 1 + j) + a'xit — 1 + j)) + 1 + j) + a'i^{t — 1 + j)) 'j dt 

( 4(1 - X)a'i^it)a'^{t) 


= 2 


+ X {ai^{t — 1) {a\{t + 1) + a'xit — 1)) + a)j(t) {ax{t + 2) + a\{t))) j dt 

+00 1 

+ 2^ / + j)(4(l — A)Q;'^(t + j) 

+ X {2a\{t + j) + a'^{t — 2 + j) + a'^{t + 2 + j)) ^ dt. 

By Lemma 1, for t G (0,1) the values of a'/^(t — 1) and ajj(t) are determined from the initial data, 
and since they are equal to zero, we have a'/^(t — 1) = ah(t) = 0. This yields 

+00 

Lx{h) = 2^ / a'h{t + k)U{l - X)a'y^{t + k) 

+ A {2a'^{t + fc) +«(,,(! — 2 + /c) + a)^(t + 2 + fc)) ^ dt. (43) 
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If A > 0 then the roots of the characteristic polynomial p\ defined by (42) are za and In 
particular, we have p\{z\) = 0. Note that, by Lemma 1, for t G (0,1) the values of a'^{t — I) and 
of are determined from the initial data. By (35), we have 


+ 1 ) = -a'xit - 1 ) + 


(44) 


for t G (0,2). Using the representation (21) of u^{t) for t G (0,1), and using (34), we infer that 

+ = z\a\{t—l) fort G (0,1). Similarly, using (35), we have a\{t + 2) = —a\{t) + u'^{t + l) 

for t G (0,1). Using the representation (21) of u^{t) for t G (1,2), and using (34), we infer that 

a'xit + 2) = z\a'x{t) for t G (0,1). It follows that a'x{t + 1) = z\a'x(t — 1) for t G (0,2). By 

induction, using (44), (22) and (21), we get that 


(45) 


a'xit -l + 2k) = z’la'xit - 1), 

for every t G (0, 2) and every fc G N. Therefore, we have obtained that 

Xo-xit 4“ 2 T fc) -f (4 — 2X^cxx(J^ T ^) T AcrA(t — 2 -t- k'j = pxi^zx^otxit — 2 -t- fc) = 0, 

for every t G (0,1) and every /c G N*. We conclude that Lx{h) = 0. Concerning the value of the 
objective functional of (OCP)“, we infer that 

^+oo 


f + OO 

J (^{1 - X) + Xu{tfj dt 

/• + 00 , . 

> [il-X){d,yTit,0)f + XuT{tf)dt, 

with a strict inequality whenever h ^ 0. It follows that is the unique optimal solution of 
(OCP)“ for A > 0. For A = 0 the result also follows from the representation (43) of Lo{h), with 
the difference that, in this case, the characteristic polynomial is po(z) = z having the unique root 
zo = 0 . 

The inequality (23) follows from (45), since for the optimal state we have the energy 

/ (d^y^(t + 2k,x)f+ (dty^(t + 2k,x)fdx= / a'xisf ds 

Jo Jt+2k-l 

= / a'xis + 2k)^ ds = jzxi^^ / a'xisf ds, 
Jt-i Jt-i 


for every t G (0, 2) and every fc G N. Theorem 1 is proved for T = +oo. 

Remark 5. The computation of the solution with the characteristic polynomial px is related to 
techniques used for linear difference equations, or for finite-dimensional linear systems with tridi¬ 
agonal matrices (see [13]). 


4.3 Proof of Theorem 2 

We assume that T = 2n. Using (25), we have ^ and then, using (16) and the represen¬ 

tation (21) of tt“, we infer (26). We have 

{^j ^ (1 + ^a) (llj/olU^lo.i) + l|yi||L2(o.i)) ■ 
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The inequality (27) follows similarly, using (15) and (22), since 


“r-<IL2(2Mfe+2) = (/ 


Theorem 2 is proved. 


5 Conclusion 

We have discussed the influence of the objective function and of the time horizon on optimal Neu¬ 
mann boundary controls for the ID wave equation. If the objective function is the control norm 
(A = 1) and if the terminal state is prescribed exactly, then the control action is distributed uni¬ 
formly over the whole time horizon, and coincides with the control given by the Hilbert Uniqueness 
Method, which is periodic. 

In contrast, if the objective function involves an additional tracking term (A < 1), then the 
optimal control action is essentially concentrated at the starting time 0 and at the terminal time 
T, and in-between it is exponentially close to 0. We have given explicit formulas, showing that 
the control is the sum of an exponentially decreasing term and of an exponentially increasing one. 
If the time horizon is infinite (without final conditions), then only the first term remains, and the 
optimal control exponentially stabilizes the system, accordingly to the classical Riccati theory. The 
norms of the control action and of the optimal state decay then exponentially in time. These results 
show that as soon as the objective functional of the optimal control problems for the considered 
system contains a nontrivial tracking term, the optimal solution has a special behavior referred to 
as the turnpike phenomenon. 

Finally, we have shown that, if the final time T is a positive even integer, then there exists a 
weight A such that the solution of the problem of exact controllability with minimal control norm 
coincides with the solution of the infinite horizon optimal control problem along the time interval 
[0, 2]. This result justifies a receding horizon control strategy, where the first part of a finite horizon 
optimal control is used and then the procedure is updated in order to control the system over an 
infinite time horizon. 

As already said, the turnpike property has been much investigated in finite dimension (see [17] 
and references therein for a general result). In the infinite-dimensional setting, in [14] distributed 
control has been considered both for the heat equation and the wave equation. 

The turnpike phenomenon put in evidence in the present paper shows an interesting qualita¬ 
tive bifurcation of the HUM control as soon as the objective functional involves a tracking term. 
However, here, we have been able to show it by means of explicit computations. 

Several open questions are in order. 

First of all, it makes sense to consider an objective functional in which the tracking term is 
replaced with a discrepancy between the solution and a time-independent function, which is not 
necessarily a steady-state. According to the results of [17], we expect then that the turnpike 
property still holds true, and that, in large time, the optimal trajectory remains essentially close 
to the optimal steady state state, defined as the the closest steady-state to the objective. However, 
in that case, we certainly do not have explicit formulas as derived in the present paper. Moreover, 
here we only considered a functional penalizing the normal derivative at a: = 0, and then we can 
only consider a time-independent function that is a steady-state, as said in Remark 4. But if 
instead, we were considering for instance the full norm in i7^(0,1), then we could consider in the 
objective functional a term of the form ||y(t, •) — a(-)ll^; where o(-) need not be a steady-state. 
Then, what can be expected is that, in large time, the optimal trajectory remains essentially close 
to the steady-state of the form ax that is the closest possible to the target a(-). 
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For more general multi-D wave equations, the situation is open. Even if explicit computations 
can only be done in specific cases, we expect that the turnpike phenomenon is generic within 
the class of optimal control problems for controllable wave equations, and that HUM controls 
characterized by the adjoint system develop a quasi-periodic pattern, but when characterized by a 
more robust optimality cost, then, satisfy the turnpike property. 

Another open issue is the investigation of semilinear wave equations (see [3]), for which steady- 
states may play an important role. Of course, in that case, we cannot expect that the turnpike 
property hold globally, but it should also hold as well at least in some neighborhood of an optimal 
steady-state (see discussions in [17]). 
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